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Abstract 

We quantize the 1-dimensional 3-body problem with harmonic and inverse square 
pair potential by separating the Schrodinger equation following the classic work of 
Calogero, but allowing all possible self-adjoint boundary conditions for the angular and 
radial Hamiltonians. The inverse square coupling constant is taken to be 5 = 2i/(j/ — 1) 
with ^ < z/ < I and then the angular Hamiltonian is shown to admit a 2-parameter 
family of inequivalent quantizations compatible with the dihedral Dq symmetry of its 
potential term 9z^(z^ ~ 1)/ sin^ 3i;^. These are parametrized by a matrix U € U{2) sat- 
isfying aiUai = [/, and in all cases we describe the qualitative features of the angular 
eigenvalues and classify the eigenstates under the symmetry and its ^3 subgroup 
generated by the particle exchanges. The angular eigenvalue A enters the radial Hamil- 
tonian through the potential (A — \)/r'^ allowing a 1-parameter family of self-adjoint 
boundary conditions at r = if A < 1. For < A < 1 our analysis of the radial 
Schrodinger equation is consistent with previous results on the possible energy spectra, 
while for A < it shows that the energy is not bounded from below rejecting those C/'s 
admitting such eigenvalues as physically impermissible. The permissible self-adjoint an- 
gular Hamiltonians include, for example, the cases U = ibl2,±cri, which are explicitly 
solvable and are presented in detail. The choice U = —I2 reproduces Calogero's quan- 
tization, while for the choice U = ai the system is smoothly connected to the harmonic 
oscillator in the limit 1/ — > 1. 
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1 Introduction 



The Calogero model ^1 12] of identical particles on the line subject to combined inverse 
square and harmonic interaction potential is extremely popular because of its exact solvability 
and its connections to many interesting problems in physics and mathematics. (See, for 
instance, [Sj and references therein.) The Hamiltonian of the system is formally given by 

i=l i 1=2 j = l 

After separation of the centre of mass, the N = 2 case reduces to the study of the 1- 
dimensional Schrodinger operator 

Hy = -TT-X^ + T^mujV + (1.2) 



2m dy^ 2^2' 



As is widely known ||4], the spectrum of Hy cannot be bounded from below ii g < — and 
for this reason^ Calogero assumed in his work that fl' > For selecting the 'admissible 

wave functions' he imposed the criterion that the associated probability current should vanish 
at the locations where any two particles collide. This is intuitively reasonable if the inverse 
square potential is repulsive. 

Mathematically speaking, the selection of admissible wave functions is equivalent to choos- 
ing a domain on which the Hamiltonian is self-adjoint. Concerning the 1-dimensional Hamil- 
tonian Hy, it is known (see, e.g., E] or the books [3 El) that the choice of its self-adjoint 
domain is essentially unique if (? > but there exits a family of different possibilities 
parametrized by a 2 x 2 unitary matrix if (7 < In the corresponding quantizations of the 
N = 2 Calogero model the probability current does not in general vanish at the coincidence of 
the coordinates of the particles. Heuristically speaking, if < < then the particles can 
pass through each other by a tunneling effect despite the infinitely high repulsive potential 
barrier. Since this phenomenon refers to the interaction of any pairs of particles, one may 
expect it to occur also in the N particle Calogero model. 

The purpose of the present paper is to explore the inequivalent quantizations of the 
Calogero model under the assumption 

-^<9<f- (gy^O) (1.3) 

focusing on the simplest non-trivial case of three particles. We shall use separation of variables 
to define the quantizations. To explain the main point, recall that the Hamiltonian ()l.ip can 
be written as H = Hq + Hrei, where Ho belongs to the centre of mass and 

Hrel = Hr + r'^Hn (1.4) 



^The energy can still be bounded from below if g = — but this case would require a separate treatment. 
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describes the relative motion. The relative Hamiltonian Hj-ei consists of the radial operator 




2m dr'^ 



N -2 d 
2m r dr 
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(1.5) 



and the angular Hamiltonian 




N i-l 



(1.6) 



i=2 j=l 



where r is the radial coordinate on R^~^ spanned by the relative (Jacobi) coordinates of the 
particles and Aq is the standard Laplacian on S^^"^. In effect, Calogero's solution amounts 
to constructing an orthogonal basis of the Hilbert space L^(R^~^) in the factorized form 
RE,\{,'r)rj\{yL), where VL denotes the angle coordinates on the sphere 5*^"^ C R^~^ and 



This is equivalent to specifying self-adjoint domains for the Hamiltonians and ifr,A- 

It is well-known (see also ()A.7j) ) that for angular eigenvalues A < |^[3 — (A^ — 2)(A^ — 4)] 
the self-adjoint domain of H^^x is not unique [HIHlini- (For iV = 3 with the conventions ()2.4j) 
adopted later this means A < 1.) The spectra of the possible self-adjoint versions of if^.A have 
been analyzed in a recent paper ^U], which uses the eigenvalues of Hq provided by Calogero 
as input. What we do here is different, since we aim to tackle the problem of constructing 
inequivalent quantizations for Hq, as well. On the one hand, the inequivalent definitions 
of the radial Hamiltonian correspond to non-trivial contact interactions when all particles 
collide with each other (at the r = location). On the other hand, the inequivalent self- 
adjoint domains of Hq^ can be obtained by imposing different admissible boundary conditions 
for the angular wave function at the singularities of the potential term in ()1.6p . which occur 
at the collisions of any pairs of particles. This latter poses a more interesting problem to us, 
but it is also more difficult technically. We shall stick to the = 3 case for which the angular 
Schrodinger equation is still an ordinary differential equation on the circle S^. 

The content of the present paper and our main results can be outlined as follows. After 
fixing the general framework and conventions in Section 2, the self-adjoint extensions of the 
angular Hamiltonian are described in Section 3 (with some details deferred to Appendix A). It 
is explained that the most general local, self-adjoint boundary condition is given by eq. ()3.14|) . 
where Ug G U{2) determines the connection condition for the wave function at each of the six 
coincidence angles, 9 E S ()3.5j) . of two particle positions (see Figure 1). We then show that 
the dihedral Dq symmetry of the angular potential, which includes the 5*3 subgroup given 
by the permutations of the 3 particles, is maintained if the connection condition is chosen 
uniformly with Ug = U independent of 6, subject to ()3.15|) . The so-arising 2-parameter family 
of angular boundary conditions, defined by U in ()3.18|1 . has not been considered before for 
the A^ = 3 Calogero model. Calogero's study [T] corresponds, in effect, to the U = —I2 case. 

Section 4 is the heart of the paper, containing a detailed description of the eigenstates of 
the angular Hamiltonian based on their classification under the Dq symmetry. There are two 



Hnr]x = Ar^A 



Hr,xRE,x = ERe,x with Hr,x = Hr + Xr 



-2 



(1.7) 
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qualitatively different cases, according to whether U in ()3.18|) is diagonal or non-diagonal. 
In the former case the eigenf unctions over the angular sectors corresponding to the various 
orderings of the particles can be chosen independently, like in Calogero's original quantization. 
If U is non-diagonal, then the sectors are connected since the probability current associated 
with the angular Hamiltonian does not vanish in general at the coinciding particle positions. It 
turns out that the angular eigenfunctions can be written down explicitly once the eigenvalues 
are determined. The eigenvalues are shown to have the form A = (3/i)^, where /x is any real 
or purely imaginary solution of eqs. ()5.2|) - ()5.4p . These 6 equations (counting the signs) are 
associated with the 6 inequivalent irreducible representations of Dq, which are summarized in 
Appendix B. Note that Dq admits 1-dimensional and 2-dimensional irreducible representations 
and the associated eigenstates and eigenvalues are called type 1 and type 2, respectively. 

Because of the Dq symmetry of the Hamiltonian, it is consistent to regard the 3 interacting 
particles as indistinguishable. The 6 inequivalent representations of Dq become pairwise 
equivalent upon reduction to the S3 C -D3 subgroup generating the permutations of the 3 
particles. Thus the Hilbert space of the system decomposes as an orthogonal direct sum of 
3 subspaces containing, respectively, the exchange even (bosonic) states, the exchange odd 
(fermionic) states, and the states associated with the 2-dimensional representation of 5*3 (that 
corresponds to parastatistics). On the basis of cluster separability [n| or the assumption of a 
complete set of commuting observables ^21; one may truncate the Hilbert space to its bosonic 
or fermionic subspace in physical applications. We stress that our inequivalent quantizations 
are new even after such truncations. The question of statistics in Calogero models is an 
intriguing issue, as is clear from ^Hj relating the model (for u = in to fractional 

statistics. See also the review ^3] and references therein. 

The angular eigenvalue equations ()5.2|l - ()5.4j) cannot be solved explicitly for general U, 
and Section 5 is devoted to the qualitative analysis of their solutions. We provide (with some 
details contained in Appendix C) a characterization of the shape of the functions entering 
these equations, as illustrated by Figures 2 and 3, which allows us to see the general features 
of the angular eigenvalues. It turns out that, in addition to the always existing unbounded 
infinite series of positive eigenvalues, for certain boundary conditions the angular Hamiltonian 
possesses a finite number of negative eigenvalues, too. In Section 5.4 some remarks are 
given concerning the stability of the boundary conditions admitting a negative eigenvalue 
and similarly for the ones possessing purely positive spectra. 

The angular boundary conditions admitting a negative eigenvalue are physically not per- 
missible, since if A < then the spectrum of the radial Hamiltonian ()6.2|) is not bounded 
from below. This is demonstrated in Section 6, where we also describe the 1-parameter fam- 
ily of the radial boundary conditions that arise for A < 1, and characterize the corresponding 
energy spectra qualitatively. The eigenvalues of the radial Hamiltonian are found as the so- 
lutions of equation ()fi.l(i|l . where the shape of the function Fa is proved to be as illustrated 
by Figure 4. The results derived in this section are consistent with the previous analysis of 
the radial Hamiltonian ^LO^, which used a different method and adopted the (then always 
positive) angular eigenvalues from Calogero's paper. 

We can explicitly write down all eigenvalues and eigenvectors of the angular Hamiltonian 
in the four special cases U = ±l2, ±(Ti. The complete solution of the N = 3 Calogero model 
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in these cases is displayed in Section 7. As already mentioned Calogero's solution corresponds 
to t/ = — 12, but we argue that the U = ai case is in a sense more natural, since under this 
boundary condition all eigenvalues and eigenstates are smoothly connected to those of the 
standard 2-dimensional harmonic oscillator arising from Hrei ()1.4j) - ()l.fi|l in the g ^ limit. 

Some further comments on our results and on open problems are given in Section 8. 

2 Separation of variables 

Although we will consider in detail the = 3 case only, is worth explaining our strategy to 
construct a self-adjoint operator from the formal expression Hrei ()l-4p in general terms. 

The first step is to choose a domain for the formal operator so that it yields a self- 
adjoint operator of the Hilbert space L^(S'^~^), which we denote here by Hq. We assume 
that L^(5'^~^) can be decomposed into a direct sum of the eigensubspaces of Hq, 

^ ^^y^^ ^2.1) 

where A is the corresponding eigenvalue. Since L'^{'R,^~^) = (R_|_, r^^^dr) C?) L^(S'^~^), the 
direct sum ()2.1|) induces the decomposition 

L2(R^-1) = (S>x L2(R+, r^-2rfr) ® Vx. (2.2) 

The next step is to construct a self-adjoint radial Hamiltonian, Hr,x, of the Hilbert space 
L^(R+, r^~^(ir) out of the formal expression Hr^x that appears in p.7j) . Finally, we obtain a 
self-adjoint version of the formal relative Hamiltonian Hrei ()1.4|) by the infinite direct sum 

Hrei = ®xHr,x(S)idv^, (2.3) 

in correspondence with the direct sum decomposition ()2.2|1 of the Hilbert space of our system. 
Let us now specialize to A^ = 3. Following closely the lines of Calogero jlj we set 

h = 2m = l, (2.4) 

and introduce polar coordinates (r, 0) on the reduced configuration space R^ in such a way 
that we have 

Xi — X2 = r^/2sm(j) 

r 2 
Xi - Xz = r\/2 sin(0 + -vr) 

o 

r 4 

X3 — Xi = rv 2sin(0 + -vr). (2.5) 

The angle counts modulo 27r and r > 0. The angular Hamiltonian acts on functions on 

by 

^^=^^ = -;£ + |-^- (2.6) 
d(p^ 2sm''30 



The radial Hamiltonian associated with an eigenvalue A of M acts on functions on R+ by 

^ 1 3 2 2 f2 7) 

r dr 8 

Because of ()2.4j) . our assumption p.3|l on the range of g becomes 

13 , , 

-2<9<2. (2-8) 

and it will be convenient to parametrize g as 

g = 2i^{i^-l) with ^<^<| (i^T^l)- (2.9) 

In the above we have assumed that the eigenvectors of Hq form a complete set yielding 
the decomposition ()2.1|) . We shall see later that this assumption is satisfied for all self-adjoint 
angular Hamiltonians if = 3, since these operators all have pure discrete spectra. As for 
the radial Hamiltonian Hr,x, the discreteness of its spectrum follows from general theorems for 
any A^. This implies, by ()2.2j) . that the eigenvectors of Hrei ()2.3|1 obtained by the separation of 
variables span the Hilbert space L^(R^~^) whenever Hq admits a complete set of eigenvectors 
in L^{S^'^). In the subsequent sections we describe the possible operators Hq and Hr,x that 
arise for = 3, and characterize their spectra and their eigenvectors. 



3 The definition of the angular Hamiltonian M 

We now begin to study the angular Hamiltonian M given by ()2.6|) with ()2.9|) . At the formal 
level, the differential operator M admits Dq symmetry, i.e., the geometric symmetry of the 
regular hexagon. We wish to maintain this symmetry in our quantization of the Calogero 
model, and for this reason we select a self-adjoint domain for M which is left invariant 
under the Dq transformations. The physical motivation for the Dq symmetry comes from 
the following two assumptions. First, the 3 particles be identical. Second, the pairwise 
collisions of 2 particles be equivalent for the situations when the spectator particle that does 
not participate in the collision is located to the left or to the right of the point of collision on 
the line. Technically, the first assumption leads to the usual S3 symmetry group generated by 
the exchanges of the particles, while the addition of the second assumption renders the total 
symmetry group to be Dq containing the 5*3 as a subgroup. In other words, the symmetry 
generators of Dq that are not in the 'exchange- 5*3' subgroup are required in order to ensure 
the identical nature of those singular points of M at the quantum level which are not related 
by the particle exchanges, like the singular points at = and at = vr (see Figure 1). 

Referring to Figure let us note that the dihedral group Dq is generated by the reflections 
P3 and i?3 that operate on the circle as 

P3:(P^-cf), R3 : (j) ^ ^ - (j) (mod27r). (3.1) 
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Figure 1: The angular configuration space S^, with the six singular points located at S ()3.5|1 
and the six 'sectors' between the consecutive singularities (left), and with the axes of the 
reflection symmetries of the angular Hamiltonian (right). 



Their composition is rotation by |, 

7^z =i?3oA:0^0 + ^, (3.2) 
which generates the cyclic subgroup Cq C Dq. The other reflection elements are provided by 

p„ = (7^|)"oP3o(7^|)-^ i^„ = (7^|)"o/^3o(7^z)-" for n = i,2. (3.3) 

The reflections Pi generate an 5*3 subgroup of Dq, which we call the 'exchange-5'3', since it acts 
by permuting the original particle positions. For example, P3 exchanges Xi and X2 fl2.5|) . The 
reflections Ri generate another S3 subgroup, which we call the 'mirror-S'3' in what follows. 
The intersection of these two 5*3 subgroups is generated by the cyclic permutation, V, of the 
particles given by 

V := (7^|)-2. (3.4) 

The Dq transformations map to itself the set of angles, S, corresponding to the singular points 
of M 

5 = {—I A; = 0,1,2,3,4,5}. (3.5) 

For any function ip on or on \S and g E Dq, define the function gip by 

m{<P):=i;{g-\<P)). (3.6) 

Obviously, ^ is a unitary operator on the angular Hilbert space L'^{S^). It will yield a sym- 
metry if it maps the chosen self-adjoint domain of M to itself. 

We apply a fairly well-known procedure (described in great detail for example in j^) to 
specify self-adjoint domains for the angular Hamiltonian. As usual, we denote the differential 
operator M applied on some domain V by Mjy, and start by considering the 'minimal domain' 
Vq consisting of C°° complex functions on \S with compact support. The domain of the 
adjoint of Mj)^ is, in fact, the maximal domain T)i C L'^{S^) on which M 
differential operator. This means that Vi consists of those complex functions ijj on S^\S that 
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together with their first derivatives are absolutely continuous on any closed interval (in 0) 
contained in S*^ \5 and for which both ip and Mip belong to L'^{S^). It is a standard matter to 
show that = M-p^. Therefore M-p^ is a symmetric operator, and its self-adjoint extensions 
are restrictions of Mx>^ that can be obtained by imposing suitable boundary conditions on the 
wave functions at the singular points of the potential. 

Note that the deficiency indices [ZIEIE] oi Mj)^ are (12, 12) since any eigenfunction of M 
on any of the six 'sectors' on the circle (see FigurelH) is square integrable due to ()2.9|1 . One can 
check this (well-known) result by means of the explicit description of the eigenfunctions given 
in Section 4.1. The facts that Mxjq has finite deficiency indices and its self-adjoint extension 
constructed by Calogero possesses pure discrete spectrum permit one to conclude (see, e.g.. 
Theorem 8.18 in ^2]) that all self-adjoint extensions of possess pure discrete spectrum, 
and thus also a complete set of eigenvectors in L'^{S^). 

We shall describe the self-adjoint boundary conditions in terms of certain 'reference modes' 
defined in pointed neighbourhoods of the elements of S ()H.5|1 . We first choose two reference 
modes around = 0, which we denote as for i = 1,2. These are some real eigenfunctions 
of M in some neighbourhood of = (excluding = 0), normalized by the Wronskian 
condition 



After having chosen the reference modes around 0, we introduce reference modes (p^ around 
any 6 E S hj requiring that 



^(0) = (-l)V'(^i0) VA; = 1,2, ^= 1,2,3, OeS. 
This defines the iff. uniquely. We remark that 

Moreover, if the initial reference modes are chosen to satisfy 

di-<l>) = (-i)V.(0), 

then we also have 

^k''i<p) = i-^)V,{m and ^TH<i>) = d{<p-l)- 



(3. 



(3.9) 



(3.10) 



(3.11) 



Using that the reference modes are square integrable due to ()2.9j) . one can show (for 
elementary arguments, see [Hllini) that the following 'boundary vectors' are well-defined for 

any -0 G Vi. 



Bern :-- 



and 



Bern :-- 



27r 47r , , 
for 6 = 0,—,—, (3.12) 



for ^=3,^,y (3.13) 
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Here, W[iIj, v^f]e± '■= linie^±o ^lip, ^Pi]{0 + e). According to the general theory of self-adjoint 
differential operators [7|, the components of the boundary vectors give a basis of the so-called 
'boundary values' for the operator Mxjj, and the self-adjoint boundary conditions require the 
vanishing of appropriate linear combinations of the boundary values. 

We restrict ourselves to self-adjoint boundary conditions that are local in the sense that 
they do not mix the boundary values associated with different singular points of M. In fact 
(see inHmHH] and Appendix A), these boundary conditions can be written as follows: 

{Ue - WBei^p) + iiUe + WB'ei^P) = WOeS, (3.14) 

where Uq G U{2) are arbitrary unitary matrices. This local boundary condition ensures that 
the quantum mechanical probability current on remains continuous at any point of S for 
the admissible wave functions selected by (j3.14|l . 

It is important to observe that if Ue = U for all 6 & S i^. ^) and some U G U{2), then 
the boundary condition \3.14\ ) compatible with the mirror-S^ symmetry. If in addition the 
reference modes are chosen according to and the 'connection matrix' U satisfies 

U = aiUai, (3.15) 

then the boundary condition is compatible with the full Dq symmetry group. 
Indeed, the first of the above statements is a consequence of the identities 

Be{R,tP) = BrA^), B's{Ri^) = B'^^^i^), (3.16) 

which are easily verified. Under the assumption ()3.10|) . one also obtains the identities 

BeiP,^) = a.Bp^eW, B',{P,^) = a^B'p^,{^), (3.17) 

which imply the second statement. 

The most general U E U{2) subject to ()3.15|) can be written in the form 

U = e'-^e'^- = e'- ( '"".'i ) ■-= ( t ^ V (3-18) 

\^ ismp cosp J \ B A J 

The self-adjoint domains T>u C Vi of our interest are given by the boundary condition ()3.14j) 
with Ue = U in ()3.18p for all 6 E S. To denote the operator defined by applying M ()2.6|) on the 
domain Vu, we use the notation to exhibit the dependence on U. In the next section we 
will fix the reference modes and investigate the so-obtained self-adjoint angular Hamiltonians 
in detail. The operator is the one denoted by Hq in Section 2, and henceforth the 'hat' 
is generally omitted from our self-adjoint operators for brevity. 

In Section 7 the condition ()3.14|) is referred to as the 'Dirichlet' case if Ug = — 12 for 
all 9 E S, the 'Neumann' case if Ug = I2, and the 'free' case if Ug = o^. To explain this 
terminology jH], note that the boundary condition of the form ()3.14|) can also be considered 
for sufficiently regular potentials for which the reference modes are smooth at Q and can be 
chosen to satisfy V9^(6'±0) = —(5^,2 and ^(6'±0) = In such circumstances ()3.14p becomes 
the standard Dirichlet condition ■?/'(^) = if f/^ = — 12, the Neumann condition = if 

Ug = I2, and the free boundary condition requiring the continuity of ^ and ip' at 9 if Ug = ai. 
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4 Eigenvalue equations and eigenvectors of M 



In this section we develop a convenient formalism to study the eigenvalue-eigenvector equa- 
tions for the angular Hamiltonian given by ()2.fij) . ()2.9j) with the boundary condition ()H.14j) 
specified hy Ue = U from ()3.18j) for all 6* G 5 The cases of diagonal ('separating') and 

non-diagonal ('non-separating') U are considered separately. In both cases the eigenvectors 
are obtained almost automatically once the eigenvalues are determined, and they are classified 
according to the representations of Dq symmetry group ()3.6|) . The eigenvalue equations will 
be further studied later. 



4.1 Preparations 

For any complex number /i, consider the functions 



t'i,/.(0) := I sm30| F \ , ' + 2' ' 

/ l - z/-/il-z/ + /i 3 .2, 



V2,,{<i>) := I sin30r-'^F -z/ + -; sin^ 30j , (4.1) 

where -F(a, 6, c; 2;) is the standard hypergeometric function. Locally, the f are linearly in- 
dependent^ eigenfunctions of the differential operator M ()2.()j) with eigenvalue [1 . They 
are square integrable in neighbourhoods of the singularities ()3.5|) . since v satisfies fl2.9p . Un- 
fortunately, these functions are not differentiable at those values of (p for which sin^ 30 = 1, 
where there is no singularity of the potential in ()2.6p . and also do not satisfy the boundary 
condition in general. Essentially, our problem is to select those linear combinations of the 
Vi^fj, that are smooth on S^\S and satisfy ()3.14j) . Since M with ()3.14|1 is self-adjoint, we can 
restrict our attention to real or purely imaginary values of /x, for which is real. Since ±/i 
define the same eigenfunctions, we can take /x to be either non-negative or of the form i|/x|. 

Below, we shall need the limiting values 



ai(/i) := lim fi,;,(0), &i(/i) := lim d^Vi^^{(f)). (4.2) 



Explicitly, we have 



(4.3) 



r(^)r(^) ' ^"^^ r(^^^^)r(^^^^; 



Now we fix the boundary condition by choosing the reference modes to be 

^?(0) = (3(2z/ - l))-^i;i,^„(0)[e(0) - e(-0)], 
^^{(1)) = -3-H2U - (4.4) 



^The independence requires ^) ^ Z, which explains why g = — ^ is excluded from our treatment H2.9|l . 
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where (j) may vary as— | < (J) ^ < |, Bis the usual step function, and /xq is an arbitrary 
real number. (The value of Ho does not affect the boundary condition; only the asymptotic 
behaviour of the reference modes around matters.) 

From this point on we study the operator := Mjy^, where the self-adjoint domain 
is specified by ()3.14j) using a matrix U from ()3.18|) and the above reference modes. Note that 
eq. ()3.1()|1 holds, and thus admits the Dq symmetry. 

It is clear that the eigenfunctions of yield smooth functions on \ S. In order to 
find them, we adopt the following strategy. We first write down all eigenfunctions of the 
differential operator M that are smooth on \ S. We then select the admissible eigenvalues 
and eigenfunctions by imposing the 'connection conditions' ()3.14j) . 

Let us define the functions !]]_ on \ S by 



and 







ri-,M) = { -a2(/i)t^i,^(f -4>) + ai(/i)^2,M(f 




if < 
if f < 
otherwise 



< ^ mod 271 



< I mod 271 



(4.5) 



if < < f mod 271 
if f < < f mod 2n 
otherwise 



(4.6) 



The functions ri]_^^ are supported on 'sector 1' on the circle (see Figure H]) , and enjoy the 
symmetry property Rsrj], ^ = itr/^ ^. In correspondence with the other five sectors on 5*^, we 
introduce the rotated functions r]^^^ by 



vlM=viM-ik-l)^), for k = 2,...,6. 



All these functions belong to C°°{S^ \ S) and are eigenfunctions of M, 



(4.7) 



(4. 



They are square integrable for any fi. The most general smooth eigenfunction of M on S*^ \iS 
can be written as the linear combination 



(4.9) 



k=l 



with arbitrary complex numbers C^. 

Our problem is to select those linear combinations ()4.9|) that are compatible with the 
boundary condition ()3.14j) . To do this, we need to compute the boundary vectors Beij]^) and 
Bg{rj^) ioT 6 E S fl3.5|) . By using the standard formulae for the asymptotic behaviour of the 
hypergeometric function, we easily find that 



B^^^ = (3(2u - 1))- 



-clbM-cUM 

C^&2(/i) + Cla2(/i) " 
Clb2{fi) - C!a2(/i) 



(4.10) 
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The boundary vectors at any 6 E S are then obtained from the identities 

■3- -J d 

which follow from ()3.16|) and ()3.17|) . For example, for rj^ in ()4.9|] 

6 

r]^on^ = Y^ ('^ivl,t, + Ctrit^^ with A; = A; + 1 mod 6 



k=l 



and thus 



i?^M = (3(2z/-l))- 



i?;(7^^) = (3(2z/-l))^ 
All other boundary vectors can be calculated similarly. 



-C|6i(^) + Clai(/i) 
Clb2ii2) + C^_a2ifi) 



(4.11) 



(4.12) 



(4.13) 



4.2 The separating cases 

For reasons that will become clear shortly, we call 'separating' the cases for which U is 
diagonal: 

■ e'" 
e^" 



U 



(4.14) 



To start, notice that the constants appear only in the upper component of the boundary 
condition ()3.14j) for = and in the lower component for 6* = By adding and subtracting 
these two equations we obtain 



[i(e'" + l)a2(^) - (e'" - l)ai(;u)]Ci = 0, 
[i(e- + l)bM - (e- - lMfi)]Cl = 0. 



(4.15) 



The equations for decouple for different values of k, and they are the same for any k. 
Let us now look for a special solution for which = for k 1 and 

|c;|2 + |ci|Vo. 

Assume that 



If both and Cl were non-vanishing, then we could conclude from ()4.15|) that 

02 (a*) ^2 (a*)' 



(4.16) 
(4.17) 

(4.18) 
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but it is possible to check (see equation ()4.34|) below) that this can never happen. Hence we 
have the following two sets of solutions: 

V, = V-,,, (4.19) 

V, = Vl,,- (4.20) 

The nature of the solutions of these 'eigenvalue equations' will be analyzed later. It is clear 
that the general solution corresponding to an eigenvalue determined by ()4.19p and by ()4.20j) 
is given respectively by 

6 6 

< = E^"^-.M and V^ = J2cUi,, (4.21) 

k=l k=l 

with arbitrary coefficients C^. 

If e'" = —1, then we obtain the following two sets of solutions: 

case A: ai{fi) = and case B: &i(/u) = (4.22) 

with the associated eiqenfunctions r]^ and of the same form as in ()4.21|) . respectively. In 
these cases the relevant formulae of t]^ ^ simplify, since in (14. 5j) and ()4.6j) only the contributions 
of vi^^ survive. The eigenvalues are given according to 

ai(/i) = ^ fi^ = {2n + l + uy, ri = 0,1,2,... (4.23) 

bi{fi) = ^ i? = {2n + uf, n = 0,l,2,... (4.24) 
This coincides with Calogero's solution pi of the angular Schrodinger equation for iV = 3. 
If e^" = 1, then the following two sets of solutions result: 

case A: 02 (yu) = and case B: fc2(/^) = (4.25) 

with the corresponding eigenfunctions provided by 1)4.211) . Now only the contributions of ^2,/^ 
survive in //^ ^ , and the eigenvalues are furnished by 

a2(/i) = ^ ^? = {2n + l + {l-v)f, n = 0,l,2,... (4.26) 

fe2(/i) = ^ /i2 = (2n+(l-z/))2, n = 0,l,2,... (4.27) 

The boundary conditions with diagonal U are 'separating' in the sense that they admit 
a basis of eigenfunctions such that each eigenfunction is supported on a single sector on the 
circle between two singular points of M ()2.(i)l . In the separating case the multiplicity of each 
eigenvalue is 6, in correspondence with the 6 arbitrary coefficients in the eigenfunctions r]^ 
and t]^ above. Let us denote by and the characters of the representations of Dq 



case A: — -— = cot — 

2 

case B: -r-7-T = cot — 

&2(/i) 2 
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defined by ()3.6|) on the eigensubspaces of M^, spanned by {rj'l^^} and {^7+^^}, respectively. It 
is straightforward to calculate that 

X^{Rk) = -xtiRk) = 2, xi{Pk)=xi{n\) = for A; = 1,2, 3, X = A,B. (4.28) 

In terms of the irreducible characters described in Appendix B, one obtains 

Xt = X'^ + X'- + X^'^ + X^^\ = + X^- + X^'^ + x}^^- (4.29) 

This fixes the decomposition of the eigensubspaces in the separating case into irreducible 
representations of Dq. If desired, one could easily implement the decomposition explicitly. 
The states associated with the characters ^^e 'bosonic' and those associated with x^~ 
are 'fermionic' under the exchange- 6*3 subgroup of Dq, for q = ±. The explicit form of these 
bosonic and fermionic states is the same that appears in eqs. ()4.54|) - ()4.57|) below. 



4.3 The non-separating cases 



Let us assume that U in ()3.18|1 is non-diagonal. Then we can rewrite the boundary condition 
for 7]^ ()4.9p in the form 



cx 



r(/i) 



cx 



with k = k + 1 mod 6, V/c = 1, . . . , 6. 



(4.30) 



The 'transport matrix' T(/i) is independent of k because of the symmetry of the problem. 
To find r(/i), one may consider ()3.14p for ip = rj^ with 6 = Q and rewrite this equation as 



where 



{h - bA) {a - aA) 
-IB -aB 



CI 



CI 



N_ 



bB -aB 

'-b + bA) {a - aA) 



and we have introduced the complex quantities 

a := ai + ia2, b := bi + 162- 



(4.31) 



(4.32) 



(4.33) 



These, and hence N±, depend on fj,, but we dropped this from the notation. The afc(/u) and 
6fc(/i) in ()4.3j) are real, since n is real or purely imaginary, and a and b are the complex 
conjugates of a and b. With the help of the refiection formula of the F-function, we find the 
useful relations 



ai(/i)62(Ai) - &i(/^)a.2(Ai) = S - Qu, ai(/i)62(/i) + 6i(yu)a2(/^) = (3 - 6z/) 
The first relation implies that 

det = det N_ = {ab - ba)B = -2i(3 - 6u)B, 



cosvr/i 



cosvrz/ 



(4.34) 



(4.35) 
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which is non-zero precisely if U is non-diagonal. Thus we have 



T(/i) = N-\fi)N.{fi), det r(^) = 1. 



Explicitly, we obtain 



with 



T(/i) 



det A^. 



(4.36) 
(4.37) 



X = -abB"^ + (aA - a)(bA - 6), 
y = d^B^ - {aA-af, 
z = W -(hA-hf. 



(4.38) 



We see from 1)4.301) that any eigenstate is completely determined if the corresponding 
coefficients C\. are given, and hence the dimension of the eigensubspaces is now at most 2. 
The C\. and yU are of course not arbitrary, the crudest condition on them being 



\c\\ 




\c\\ 


_ c\ _ 







(4.39) 



This condition must hold since the eigenfunctions of M are (smooth) functions on \ S, 
i.e., they enjoy the periodicity ?7^(0 + 27r) = ri^{(j)). We simplify the requirement ()4.39j) by 
classifying the eigenstates according to the irreducible representations of the group Dq, which 
is possible because any operator g ()3.6|) for g & Dq commutes with M^. As summarized 
in Appendix B, the group Dq admits four 1-dimensional and two 2-dimensional irreducible 
representations. In what follows, we call an eigenstate 'type 1' if it belongs to one of the 1- 
dimensional representations, and 'type 2' if it belongs to one of the 2-dimensional irreducible 
representations. The corresponding eigenvalues will be called type 1 and type 2 as well. To 
describe the states in terms of these representations, it is convenient to diagonalize, together 
with , the operator T that represents the rotation^ 

T := {n^y\ (4.40) 
The action of T on a general eigenfunction ()4.9|) is given by ()4.12|) : 

if V^l = ^ C^vi,,,, then frj^ = ^^ C^t]^^^ with k = k + 1 mod 6. 



(4.41) 



k,± 



k,± 



Since = id, the possible eigenvalues of T are the sixth roots of unity, , = 0, 1, . . . , 5, 
which are ±1, ±j, ±j in terms of the cubic root j := e~ . By combining ()4.4ip with ()4.30|) . 
we see that a joint eigenstate rj^ of and T, for which 



fri^ = Tri^, r e {±1, ±j, ±j}. 



(4.42) 



^Our subsequent treatment was inspired by |10) . where the cyclic permutation operator V = rather 
than T was used to solve the system for ^ < v <1 with U = ai, in effect. 
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is equivalently given by its coefficients that satisfy 



T 



(4.43) 



This imphes (14 .3911 but is more useful to characterize the eigenstates, because one can tell 
the 'type' of the eigenstate from the eigenvalue r: it is type 1 if r = ±1 and type 2 if 
r G { ±j, ±j} (see Appendix B). 

The existence of a non-zero eigenvector in ()4.43|1 is equivalent to 



= (detiV+)Met(r(/i) - rla) = [x(/i) - r det iV+]^ - ?/(/i)z(/i). 

Combining this with detT(yu) = 1, and using ()4.35j) and ()4.38j) we find 

x(/i) _ {aA - a){hA -b) - obB'^ _ 
detAr+ ~ {ah - ha)B ~ 

where 3?(r) denotes the real part of r. For r in ()4.42|1 



(4.44) 



(4.45) 



3?(r) 



2r 



±1 if r 
±\ ifr 



±1, 



(4.46) 



The admissible values of yU are determined by the spectral condition ()4.45|) once r is specified. 

Let us now focus on the type 1 eigenstates for which r = ±1. In this case the spectral 
condition ()4.45p can be factorized as 



(aA-a- 3?(r)ai3) {hA-h + ^{T)bB^ = 0. 



(4.47) 



This allows us to classify the solutions according to the possibilities as to which of the two 
factors in ()4.47|) vanishes and which of the signs 3fJ(r) = ±1 is chosen. These four possibilities 
are listed as 

cases A±: a{n) = {A±B) a(/i), (4.48) 

cases B±: b{fi) = {A±B) b{fi). (4.49) 

Note that we have assigned 3?(r) = =i=l to A-t and 3?(r) = ±1 to B-t, respectively. The 
conditions ()4.48|) . ()4.49p can be spelled out in the more explicit form 



cases A-I-: 
cases B-I-: 



Qi(/i) 
a2(/i) 



cot 



cot 



(4.50) 



(4.51) 



which should be compared with ()4.19p and ()4.20|) . We shall analyse the roots of the eigenvalue 
equations ()4.50|) . ()4.51|) later, and here we just mention the special case in which cot 
becomes divergent. This happens precisely if {A ± B) = e'^"^^-* = 1, whereby the above 
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equations formally give 02 (/i) = and 62 (yu) = 0, respectively. These conditions result also 
directly from ()4.48|) . ()4.49p . because they simplify as 



a(/i) = a(yu) and b{fi) = b{fi) if 



,i{"±/3) 



(4.52) 



which can be solved trivially. Similarly, we have 

a{fi) = -a(/i) and b{fi) = -b{fi) if e'^'"^'^) = -1, 



(4.53) 



and then the solution is obtained from ai(/i) = and 61 (/i) = according to ()4.23|1 and ()4.24|1 . 

Since the cases A± involve only a(/i) and a(/i) while B± involve only b{fi) and b{fi), each of 
these conditions occurs if the eigenstate r]^ ()4.9|) consists exclusively of r]^^^ or 77^^^ (see ()4.5|) 
and ()4.fij) ). By using this observation and the assignment of 3?(r) = =i=l to A± and 3?(r) 
to B-t, we immediately find the corresponding eigenstates to be 



±1 



k+l^k 



k=l 



k=l 



(4.54) 



= (4-55) 

k=l k=l 

One can check that the reflections Rk,Pk G -De act on these states according to 



PkV^^ = ±r]'^^ for X = A,B, 



(4.56) 



which fixes their association with the four 1-dimensional representations of D^. If we label 
the type 1 eigenstates by the parities analogously to the labeling of the type 1 characters of 
Dq in Appendix B, then we can write 



PL 5 



PL 5 



(4.57) 



In particular, the states r]^+ are 'bosonic' and the r]^~ are 'fermionic' under the exchange-^s. 

The eigenstates ()4.54|) . ()4.55|) may also be confirmed by examining the eigenvectors in 
()4.43j) . In fact, the conditions ()4.48j) and ()4.49|1 are equivalent to ?/(/i) = and z{n) = 0, 
respectively, and hence the 'transport matrix' becomes triangular for the type 1 eigenvalues. 



cases A-I-: T(/x) 



det A''4. 







(4.58) 



cases B 



T(/i) 



±1 



dct Af+ 

±1 



(4.59) 



The obvious eigenvector of T(/i) gives rise to the corresponding eigenstate in ()4.54|) . ()4.55|) . 
It is also possible to show that z{fi) and y{fi) can never vanish simultaneously, and thus T(/i) 
above is truly triangular. This implies that the multiplicity of each type 1 eigenvalue is 1. 
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Next, we turn to the type 2 eigenstates that span the 2-dimensional eigensubspaces of 
. The states belonging to the Dq representation of character x^'^^ (see Appendix B) may 

(2) — 

be denoted as rjl^r with r = — j, —j and those belonging to the representation of character 



v(2) 



X 



as rjli^T with r = j, j. For each 3?(r) = ±i the admissible values of fx are determined by 



the spectral condition ()4.45|1 . which can be expanded straightforwardly as 

cos /3 + cos a 02 &2 
sin P a2&i — ai&2 



sin a 0162 + ct2&i _|_ cos (3 — cos a aibi 



sin 0162 — ct2^i sin /3 0261 — ai&2 

Upon substituting ()4.34|) into ()4.60p . we obtain 



+ 



3ft(r). (4.60) 



sinacosTT/i cos /3 — cos a cos /3 + cos a 

+ T7^ ^«l(Ai)Ol(/^) + 77 nO'2\Nb2\H = ^i(t) 



sin/?cos7rz/ ' (6z/ — 3) sin/3 ^^'"^ ^^'^^ ' (6i^ — 3) sin/5 ^^'"^ ^^'"^ ^ ^ ^ 

The products aibi and 0262 do not simplify to trigonometric functions. The best we can do 
is to rewrite them using the identity 



T{z)T{z+-) 



r(2.)r(l)2^-2^ 



as 



ai{fi)bi{fi) 



6r2(z/ + |)22(--i) 

r(z. + ^)r(z/-/i) 



a2(/i)&2(/i) 



-2v 



V(l-u + ii)V{l-u- fj)' 



(4.62) 



(4.63) 



In contrast to the type 1 case, now some further work is needed to find the eigenvectors 
in ()4.43|) for any admissible /i solving 1)4.611) . However, since there must be two independent 
eigenvectors (with eigenvalues r and r) to form a 2-dimensional representation of Dq, the 
condition ()4.39|) must actually hold automatically for arbitrary C\. This implies that we 
have the matrix identity T^(/i) = I2 for any type 2 eigenvalue (as can also be confirmed by 
a direct computation). Thus, we can associate with any eigenvalue r of T[ii) the projection 
operator 

k=l 

which satisfies T(/i)7rT-(/i) = r7rT-(/i), 7r^(/i) = TTrifJ') by virtue of rf 
in ()4.43|) is provided by 

(ji ] 



(4.64) 

1. Then the eigenvector 



7 
5 



(4.65) 



with arbitrary complex numbers 7, 5 for which the eigenvector is non-vanishing. The projec- 
tion operators may be evaluated explicitly as 



7r^(/i) 



1 
2 



3 det N. 



(4.66) 



where 9'(t) G {±^} denotes the imaginary part of r. In principle, the construction is 
completed by using ()4.30)1 and ()4.9|1 to generate the type 2 eigenfunctions of . 

The eigenvalues of cannot be presented explicitly in general. In the next section we 
investigate some features of the eigenvalues for general U (I3.18jl . We know of four explicitly 
solvable cases, corresponding to U = ±l2, ±cri. These are discussed in Section 7. 
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5 Characterization of the eigenvalues of 

We have seen that the eigenvalues of the angular Hamiltonian are given by 

A = (3^)2, (5.1) 
where /x is a solution of one of the following 'eigenvalue equations'. First, 



r(^ ±|±^)r(^±|^) _ r(z/ + |) a ±(3 
^)r(^^) ~ T{-i^ + l 

r(^)r(^) _ r(z/ + i) ^ , a±/3 



Fa(u) := \ \ \ ^ = ^- tan (5.2) 

^vw r(^^)r(^^) r(-z/ + |) 2 ^ ^ 

or 



Fb(u) := — / f \^ = ^- tan -. (5.3) 

r(i^)r(i^) r(-i/ + f) 2 ^ ^ 



Second, 



sin a cosnfi cos /? — cos a cos /9 + cos a 1 

smpcosTTZ/ (oz/ — 3)sinp (oz/ — 3)sinp 2 

Equations ()5.2|1 and ()5.3|1 control all eigenvalues for the separating boundary conditions, 
for which (3 = modulo vr, and the 'type 1' eigenvalues for the non-separating boundary 
conditions. In the latter case, the states corresponding to the solutions of ()5.2|) are odd 
while those corresponding to ()5.3|) are even with respect to the mirror-Ss symmetry. Their 
parity under the exchange-S'3 symmetry is given by the ± in the argument of the tangent 
function on the right hand side. Equation ()5.4|1 governs the 'type 2' eigenvalues for the non- 
separating boundary conditions. In each of the cases, we are interested in the values of fi that 
satisfy either /i = or /i = i|/i| in correspondence with the non-negative and the negative 
eigenvalues of . In fact, one of the main questions is whether negative eigenvalues exist 
or not. 

In our subsequent analysis of equations ()5.2|) and ()5.3|) we assume that I tan ^1 < 00. If 
tan ^ diver ges or vanishes, then all solutions of the corresponding equations are real and 
can be written down explicitly, as discussed in Sections 4 and 7. 



5.1 Type 1 negative eigenvalues 

We now study the possibility of negative eigenvalues arising from equations ()5.2|1 . ()5.3j) . which 
govern the type 1 eigenvalues for non-separating U and all eigenvalues for separating U. By 
setting fj, = ix with x > 0, we consider the real functions on [0, 00) defined by 

I -p l+i^+ix ^ 1 2 |-p^ t/+ix ^|2 

X 1-^ FA(ix) = - — ^, . , ,„ and x 1— >• Fniix) = - — --; — ^, . , ,„ . (5.5) 

' j 2—v+\x j2 ^ ' jp^ 1— j2 ^ ' 

We show below that these functions are strictly monotonically increasing on [0, 00) and they 
grow without any bound as x tends to +00. Therefore there exists at most 1 negative eigenvalue 
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of that arises as a solution of ()5.2j) . for any fixed sign ± on the right hand side, and such 
eigenvalue occurs precisely if the parameters of satisfy 

Fa{0) = < tan (5.6) 

Similarly, for any fixed sign ± on the right hand side of ()5.3|1 . we obtain at most 1 negative 
eigenvalue, which occurs if and only if 

F m-^^<i^^^^tan"^^ r57) 
Fb{0) - < tan (5.7) 

To prove the above claims regarding FaIix), we inspect its logarithmic derivative: 
2 (ilogF^(ix) + + f I + u - ix\ f 2 - u - ix\ f2-v + ix 

(5.8) 

Viz) 

m = Yji- (5-9) 



i dx \ 2 

with the standard notation 



Recall the identity (|2I|: 8.363 4.), for any real ^ and r], 



^(e + i^) - ^(e - = E + (5-10) 

By using this and adding the two series, we obtain 

c;iogF^(ix) _^ 2a:(2z/-l)(3 + 4fc) 

^ [x2 + (2- z/ + 2A;)2][x2 + (! + !/ + 2A;)2]' ^ 

which is positive for x > as (2z/ — 1) > 0. Since Fa^ix) > 0, this imphes that x \—>- F^(ia;) is 
a strictly increasing function on [0, oo). Next we recall (120]: 8.328) that 

lim \T{^ + ir])\e^^i^\T]\^-^ = V2^. (5.12) 

|r?|^oo 

It follows from (jKT^ and that 

lim FA^ix) = lim ( — ) = oo, (5.13) 
as had been claimed. The analogous properties of Fsiix) can be verified in the same manner. 
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Figure 2: Fa (left) and Fb (right), as the function of /x [see ()5.2|) - ()5.3p ]. for /i > 0, with 
u = 2/3. The dashed hues are located at /i^ (jOIj) for Fa and fl5T9|l for Fb- 



5.2 Type 1 non-negative eigenvalues 

Next we describe the shape of the functions /i i-^ Fa{ij) and /i ^— Fb{ij) for /i > 0, as illus- 
trated by Figure El This permits us to see the main features of the non-negative eigenvalues 
dnH} of furnished by the solutions of and 

The function Fa fl5.2|) is smooth for yU > 0, except at the values 

/i^ = (z/ + 1) + 2m, m = 0,1,2,..., (5.14) 

where it becomes infinite. More precisely, as is easily checked using the properties of the 
F-function, Fa{ij) approaches +oo and — oo as approaches /x^ from above and from below, 
respectively. It is positive at /i = 0, and it takes the value zero at 

/xli, = (2 - z/) + 2m, m = 0,1,2,... (5.15) 

Note that 

/^m < /^m < /^m+l < /^m+l, m = 0, 1, 2, . . . (5.16) 

It can be shown (see Appendix C) that Fa is strictly monotonically decreasing for /i G [0, /i^) 
as well as for jj, G {fi'^ , fi'^^i) for any m > 0. 

Since the constant on the right hand side of ()5.2j) is finite, we conclude from the above 
that ()5.2|) admits a unique solution in the range (/im'/^m+i) (!5.14|) for any m > 0. There is 
an additional solution in the range [0,/i^) if 

Let us sketch the analogous description of the function Fb{^) for /i > 0. As is readily 
verified, Fb{ij) changes sign from positive to negative as /i passes through the zero locations 
given by 

/i° = |l-z/|, f^ = {l-u) + 2m, m = l,2,.... (5.18) 
It becomes +oo and — oo as it approaches 

/i^ = z/ + 2m, m = 0,l,2,... (5.19) 
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from above and from below, respectively. We have 



fil < fi^ < frn+i < m = 0,1,2,.... (5.20) 

Similarly to the case of Fa, one can prove that Fb is strictly monotonically decreasing for 
11 G [0, /i^) as well as for /i G (/ij^, J^m+i) /o'^ ^^?/ > 0. 

As a consequence of the shape of Fb, for any finite constant on the right hand side, ()5.3j) 
has a unique solution in the range (/i^, /i^_,_]^) ()5.19j) for any m > 0. There is an additional 
solution in the range [0,/i^) if 



5.3 Type 2 eigenvalues 



The formula of the function F2 in ()5.4p is rather complicated for general U , for this reason 
we shall be content with some remarks on the generic properties of the type 2 eigenvalues. 

As an illustration, let us first investigate the equation of type 2 non-positive eigenvalues, 

F^iix) = ±i, x> 0, (5.22) 
in the special case a := — |. Under this assumption F2{}x) simplifies as 

/• \ ^ f coshTTo; cos/3 ,7 7 w \\ r^ox 

^2 IX = r + ^ [aih + a2&2 IX . 5.23 

sm/5 \| cos7rz/| (6z/ — 3) / 

We see from ()4.63|) that afc(ix)6fc(ix) > for /c = 1, 2, and therefore 

F2{ix) > ——r^ r > 1 if < P <^ = -a. (5.24) 

sm/7| cos7rz/| 2 

Thus does not admit non-positive type 2 eigenvalues for these choices of the parameters, 
which include the explicitly solvable case /3 = —a = | for which U = ai. On the contrary, if 
> I = —a and P is near enough to |, then we have 

^2(0) > 1 and lim F2(ix) = -00, (5.25) 

X— >oo 

which implies the existence of at least 2 negative eigenvalues of for the corresponding U 
()3.18p . The second relation in ()5.25p holds since for large x the term a2{ix)b2{ix) dominates 
F2, unless it is multiplied by zero, as is easily seen from ()5.12|) . 

The above example shows that, like in the type 1 case, the existence or non-existence of 
type 2 negative eigenvalues depends on the choice of the parameters a, j3. One can check that 
limj.^oo |-^2(ix)| = 00 always holds, and hence any choice leading to 1-^2(0)1 < | guarantees 
the existence of such eigenvalues. It is also clear that there can be only finitely many negative 
eigenvalues, the maximal number we found in numerical examples is four. 
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Figure 3: F2{fj,), for imaginary values of fi = ix (left) and for real values of /i (right). The 
dashed lines lie at ±i. For the left figure, u = a = ^tt and /3 = Jjvr, a case when the 
number of solutions of ()5.4|1 is four. For the right figure, the parameters z/ = a = ^vr, 
and P = Y^'^ are chosen so that the plot exhibits the behaviour of each of the three terms 
of F2 in ()5.4|) . Initially, the rapidly decreasing middle term dominates, then the constant 
amplitude first term rules, and for larger values the third term gradually becomes the most 
significant one. Each term oscillates with the same frequency. 

Since the multiplicity of the 2-dimensional representations of Dq in L'^{S^) is obviously 
infinite, has infinitely many type 2 positive eigenvalues. We can understand the approx- 
imate behaviour of the solutions of ()5.4|1 for large positive n by using the relation (for any 
real y) 

lim T{y - fi)T{y + fi){y + nY'^y sin 7C (fx + 1 - y) = vr. (5.26) 

Since > |, this implies by ()4.fj3j) that ai(/i)6i(/i) tends to zero as /i — ^ 00, and a2(/i)&2(At) 
asymptotically equals to a constant multiple of the function 

/2(/i) := (1 -z/ + /i) 2'^" Sin 7r(^ + z/). (5.27) 

If (cos a + cos/?) 7^ 0, we thus obtain from ()5.4j) that F2{fj,) can be approximated by a 
constant multiple of /2(/i) as /x — >• cx). Therefore, for large enough /x, there exist two solutions 
of \F2{fi) \ = I between any two consecutive extrema of the function /2(/x), and actually these 
solutions lie near to the zeroes of f2- If (cos a + cos (3) = 0, then the first term of F2 ()5.4j) 
is the dominant one. If cos a = cos (3 = 0, then only this term remains, and we can find 
the eigenvalues explicitly as discussed in Section 7. If desired, one could derive more precise 
asymptotic estimates for the solutions of ()5.4|1 by developing the above arguments further. 
The behaviour of F2 near to = can be rather complicated in general (see Figure Ej) . 

5.4 On the lower-boundedness of the energy as a condition on U 

We have seen that admits, in addition to its infinitely many non-negative eigenvalues, 
a finite number of negative eigenvalues, too, for certain values of the parameters a, (3 of 
U (I3.18|l . This result is important, since — as we shall demonstrate in Section 6 — the 
existence of a negative eigenvalue of implies that the energy spectrum of the model (i.e.. 
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the spectrum of Hrei (ESI)) is not bounded from below. These cases are to be excluded 
in hypothetical physical applications. It is complicated to precisely control the conditions 
for the non-existence of negative eigenvalues of M^, but there are certainly cases for which 
this occurs. For example, if a = — | and < /? < |, then tan^y^ < 0, and thus neither 
the inequalities ()5.fi|l . ()5.7|1 nor the corresponding equalities can be satisfied. Together with 
()5.24j) this implies that the spectrum of is positive for these choices of the parameters of 
U fl3.18|) . The same holds, for instance, for the explicitly solvable cases U = ±l2, ±cti. 

Next, we present a stability result concerning the positivity of the spectrum of . Sup- 
pose that all eigenvalues of are positive for some U = f/(ao,/3o) and the parameters 
(ao, Po) are generic in the sense that they satisfy the following inequalities: 

I tan ° ^ I < oo, sin /?o 7^0, (cosao + cos/5o) 7^ 0. (5.28) 

Then it can be proven that with U{a,l3) has only positive eigenvalues for any 

parameters (a,/5) near enough to (aoi/^o)- 

To verify the above statement, we first observe that neither any of the inequalities in 
()5.6|) . ()5.7|) nor the corresponding equalities can hold for (a,/5) near to (q;o,/?o) by continu- 
ity. Therefore all type 1 eigenvalues must be positive for such parameters. To exclude the 
possibility of type 2 non-positive eigenvalues, which would arise from the solutions of ()5.22|1 . 
notice that F2 fl5.4|l can be written in the form 

F2{ix,a,(3) = e"^ [Ko{a, f3)Ko{x) + K.ia, (3)x^-^''K_{x) + K+{a, f])x'^''-^K+{x)] , (5.29) 

where the functions Ka satisfy 

Kaix) > Vx > 0, lim Kaix) = 1 (Va G {0, ±}). (5.30) 



We see from ()4.63p and ()5.12j) that the above relations are valid with 
and 

1 sin a 



Ko(«,/?) 



2 cos vrz/ sin P 



^ r2(i. + l)22(-i)cos/3-cosa 

K-{a,/3) = — , 

7r(2z/ — 1) sm p 

r^(f-z/)2-^-cos/3 + cosa 

K+[a,(3) = — . 5.32 

7r(2z/ — 1) sm p 

Let us assume that ()5.28|) is satisfied and n+{ao, /?o) > (the case of the other sign is similar). 
We can choose a neighbourhood Vi of (ao, Po) and constants 7a > so that 

K+(a,/3)>7+, |fi:o(a,/3)| < 7o, |K_(a, /3)| < 7_, \/{a,P)eVi. (5.33) 
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Then we fix some Xq > for which 

F2{ix,a,(3) > e"^ [-f+x'^''-^K+{x) - ^oKo{x) - -f^x^-^" K^{x)] > 1, Va; > xq (5.34) 

and W{a,P) G Vi. This ensures that there is no solution of ()5.22p for x > xq. Since Kq{x) 
and x^'-^'^-^^Kiix) are bounded on [0,a;o], for any e > we can find a neighbourhood V2 of 
(ao, Po) so that 

|F2(ix,a,/3) -F2(ix,ao,/3o)| < e if {a,p)eV2, < x < xq. (5.35) 
By choosing e appropriately, for instance in such a way that 

F2{x,ao, Po) — -> 2e for < x < xq, (5.36) 

we conclude that ()5.22|) has no solution if (a,/?) G l^i fl V2. This implies the stability result 
that we wanted to prove. 

We can establish a counterpart of the above stability result concerning the 'impermissible' 
boundary conditions as well, for which negative eigenvalues of the angular Hamiltonian exist. 
Namely, if has one or more negative eigenvalues, then generically this property is stable 
under arbitrary small perturbations of U in ()3.18j) . Indeed, this is the case obviously if 
admits a type 1 negative eigenvalue or a type 2 negative eigenvalue which is generic in the 
sense that it arises from the graph of a; 1— *• |F2(ia;)| properly intersecting, not just touching, 
the horizontal line located at |. 



6 The radial Hamiltonian 

Recall from ()1.7|) that, at the formal level, the radial Hamiltonian reads 

^ (P l(i^322_j_'^ (6 1) 

After having characterized the qualitative features of the eigenvalues A = (3/i)^ of M^, we 
below analyze the energy levels of the relative motion of the three particle Calogero system 
defined by the eigenvalues of the possible self-adjoint versions of the radial Hamiltonian. 

Since H^^x has to be self-adjoint on a domain in L^(R+, rdr), it is more convenient to deal 
with the equivalent operator 

nr,x:=V^oHr,xo^ = -^ + liu\^ + ^^, (6.2) 
y/r dr"^ 8 

which must be self-adjoint on a corresponding domain in L^(R+, dr). It is easy to check that, 
for any eigenvalue, both of the two independent eigenfunctions of the differential operator Tir^x 
are locally square integrable around r = if and only if A < 1. For this reason |Zl El IH], 'Hr^x 
admits inequivalent choices of self-adjoint domains if and only if A < 1. It follows from general 
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theorems, collected in Appendix A from [7], that any self-adjoint version of Tir^x possesses 
pure discrete spectrum. 

If A > 1, then the unique self-adjoint domain of 'Hr,\ consists of those complex functions 
p on R_|. for which p and p' are absolutely continuous away from r = and both p and 
Ti-r^xP belong to L^(R+, dr). It is straightforward to show that the spectrum is given by the 
eigenvalues 

Em,x = 2c{2m + l + ^/X), c := W |cu, m = 0,1,2,..., (6.3) 

V o 

with the corresponding eigenfunctions 

PmAr) = r^+-^e-5-^L;f (cr^), (6.4) 

where is the (generalized) Laguerre polynomial [201; > 1. This result is contained, 
for example, in jUH]. 

From now on we consider the case 

A < 1. (6.5) 

Let ifi and (f2 be two independent real eigenfunctions of T-Cr,x associated with an arbitrary 
real eigenvalue. In fact (see [3111]), in addition to the same properties they have for A > 1, 
the functions p in a self-adjoint domain of Hr,x must now also satisfy a boundary condition of 
the following form: 

W[p,ipi]o+ 

7777 ^ = t^W, (6-6) 

where k{X) is a real number or is infinity. Here k{X) = means that v5i]o+ = 0, and 
similarly Vr[p, y)2]o+ = is required if k{X) is infinite. Our notation emphasizes that one 
can in principle choose different constants on the right hand side of fj6.6|) for different A. We 
remark that the condition ()6.6|1 can be regarded as a special case of the boundary conditions 
of the form in ()3.14|1 . where one restricts to the positive side of the singular point r = of 
Ti-r^x considered on R (accordingly U reduces to a phase), prohibiting the particle from going 
into r < from r > 0. With the 'reference modes' cpk fixed subsequently, the self-adjoint 
radial Hamiltonian specified by condition ()6.6|) is denoted as 'Hr,x,K{x)- (In the notation used 
in (j2ISl), one may substitute i^r,A,K(A) '■= ^ ° ^r,A,K(A) o Vr for Hr,x-) 

In order to determine the spectrum of Hr,x,K{x), we first write down the solutions of 

nr,xp = Ep (6.7) 

for any real number E. To do this, it is convenient to introduce 

2 E v^+1 
a := cr , ^ •= ^ ^ (6-8) 

with c given in ()6.3p . Then one can check that, if A 7^ 0, two independent^ solutions of ()6.7p 
are provided by the functions 

PE,i{r) = a5(Hv^)e-|-$(-^, ^A + 1, a). 



^To save space, we henceforth exclude the A = case from our investigation, since it would require a 
separate treatment and the final result for A = is expected to be similar to that for any < A < 1. 
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2(r) = a^(^-^)e-^'^<l>(-e - VA, 1 - TA, a), (6.9) 

where ^{a,b,z) is the confluent hypergeometric function, also called Kummer's function. 
Here, as before, we insist on the convention that either a/A = 3yU > or its imaginary part is 
positive. Up to a multiplicative factor, there is a unique linear combination of the functions 
in ()fi.9|l that lies in L^(R+, dr), given by 

r(i - Ta) r(i + Va) 

In fact, with an irrelevant factor C(V^), one has 

PEir) = C(v^)a5(Hv^)e-i-[/(-e, 1 + a), (6.11) 

where the function U satisfies U{^,1 + V^, o") = a^[l + 0(cr^^)] as a — oo (see eq. 13.1.8 
in [21]). This guarantees the square integrability of pe- For E to belong to the spectrum of 
'Hr,x,K{x), Pe niust satisfy the boundary condition ()6.6j) . 

Let us now suppose that 

< A < 1. (6.12) 
In this case the pE,k are real functions, and we fix our reference modes to be 

■= pEo,k (6.13) 
with some arbitrary real Eq. For arbitrary E and Eq, an easy calculation yields 

W[pE,k, PEo,i]o+ '■= lim (^PE,k ^^^^'^ - P^o,i-^l^^ (^) = -'iek^i'/Xy/c, (6.14) 



where e^,; is the usual alternating tensor. Therefore we find that 

W[pE,ViU T{l + VX)T{-^-VX) 

w[pE, v2]o+ r(i-VA)r(-o 



(6.15) 



By substituting back ^ fl6.8|) . we obtain the following condition that determines the eigenvalues 
of'Hr,x,K{\) under \6. VA): 

f,(,):=£Ll±i^ = _£(z^,(A) (6.16) 

Next, we analyze the shape of the function Fx, which is illustrated by Figure 0] 
Let us start by observing that F\ has zeroes at 

d = ^ + Vm = 0,1,2,... (6.17) 
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J 
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/ 









Figure 4: Fx, as the function of e [see ()6.16|) ]. for A = 1/5. The dashed hues, where the 
function diverges, are located at ()6.18|) . 



and it becomes ±oo at ej^ =F for 



+ m, Vm = 0,1,2, 



(6.18) 



By considering 



F.ie) ' 2 
with the notation ()5.9p . wee see immediately that for 

- oo < e < 



(6.19) 



(6.20) 



Fx(e) > and F'^{e) > 0. This and the asymptotics of Fx imply that Fx{e) monotonically 
increases from to +cxd as e varies from — oo to e^. If e > e^, we use the reflection formula 
!2]) for to write 



We have 
If 



^(e + ^ — ) - ^(e + ^ — , 



TT 



l-^A, , l + ^A, 
cot 7r(e H ) — cot 7r(e H J 



Wm = 0,1,2, 



(6.21) 
(6.22) 

(6.23) 

then both differences in the square brackets in 1)6.211) are negative. Since under ()6.23j) Fx{e) < 
0, it follows that Fx is monotonically increasing in this domain. If 



e°° < e < e° 



(6.24) 
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then Fx{e) > and the differences in ()6.2f|) have opposite signs. We can show that F'^{e) > 
by using the integral formulae ()C.12|) and ()C.13|) in the same way as in Appendix C, which 
proves that Fx is increasing in this domain as well. 

Supposing that k(A) ^ {0, oo}, it is clear from the shape of the function Fx that there 
exists precisely one positive eigenvalue, ^, in each interval (em;^m+i) for any non- negative 
integer m. Moreover, one obtains at most one negative eigenvalue, which occurs precisely if 

Fx{0) = ^it^ > _£iz^^(A) > 0. (6.25) 

There is also a non-neg ative eigenvalue in [0, e^) if Fa(0) < -^J^k(A). 

The eigenvalue equation (jfi.lfjj) is explicitly solvable if k(A) = or /t(A) = oo. In the former 
case the eigenvalues have the same form as in ()6.3|) . and the corresponding eigenfunction 
reduces to pe,i in fj6.9|) . which has the same form as ()6.4|) up to an irrelevant constant. If 
k(A) = oo, then we find the eigenvalues 

Em,x = 2c{2m+1 - y/X), m = 0,1,2,... (6.26) 

Under E = Em,\, Pe ()6.1()j) is proportional to pe,2, which (up to another irrelevant factor) 
gives the eigenfunction 

PmAr) = r'^~^\-'^'^"L-J~\cr^), (6.27) 
Note that all energy levels are positive in these cases. 

Our spectral condition ()6.16|) is consistent with the result in jTUj^, where the inequivalent 
quantizations of the radial Hamiltonian Hj. x fjl.7|) were considered by using a different method, 
with < A < 1 taken from the eigenvalues of two special self-adjoint versions of the angular 
Hamiltonian Hq, fll.6|l that are well- understood for any due to Calogero In our = 3 
case-study those correspond to with U = — 12 or U = I2. 

Let us now deal with the case when 

"\/A = 3/i = ix with some x > 0. (6.28) 

Then the functions pE,k are complex, and their real and imaginary parts are also solutions of 
()6.7|) . We choose our reference modes to be (up to a factor) the real and imaginary parts of 
Peo,i for some real Eq. Explicitly, with ^0 = ~ ^^2^; "we define 

tfi{r) := r^e'^" [3?<l'(— ^0, 1 + ix, a) cos(xlogr) — 53<l'(— ^0, 1 + ^x, a) sin(xlogr)] 

(p2{r) := r^e~5°" [s)^^<|)(_^g^ 1 -|- ix, a) sin(a;logr) + Q=$(— ^0, 1 + ix,a) cos(a;logr)] . (6.29) 



These eigenf unctions of Ti^.A are independent, since we find 



W[ipuV2]{r) = lim W[ipu V2]{r) = x. (6.30) 



^In WJ]I the shape of Fa in was iUustrated by a Mathematica plot, without presenting a proof of its 

properties as supphed above. 
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One also readily calculates that 



lim W[pE,i,'fi]{r) 

r— >0+ 


1 

= —ixc^ 




lim W[pE,2,fi]{r) 


-(- 


—ix) 


lim W[pE,i,<^2]{r) 

r— >0+ 


-(-H 


hix) 


lim W[pE,2,^2\ir) 

r^0+ 


-f-- 
= a;c2'-2 


-ia;) 



(6.31) 

By inserting these into ()6.6|) for p = pe ()6.1Up . we obtain the following eigenvalue equation: 

1 + ia;. 



cot arg C2'^r(l - ix)r(-e + 



-KiX) 



Ac' 



(6.32) 



Equivalently, we have to solve 



3^^ c2'^T(l - ix)r(-e + 



l+ix- 



3 C2'^T(1 -ix)r(-e + 



-k{\). 



(6.33) 



Because of the shape of the potential in ()6.2j) . one expects to find infinitely many solutions 
for e around +oo as well as around —oo. 

Indeed, we can easily prove that the accumulation points of the spectrum of 'Hrx^-W 
precisely ±oo. For this purpose, it is advantageous to rewrite ()6.32|) as 



fi(e,x) := argr(-e + ^ + 



i!}{x, k(A), c) mod 71, 



where 



X 



'i9(x, k{X), c) = arccot(— /t(A)) — arg r(l — ix) — — log c. 



In ()6.34j) we can take ^(e, x) to be the smooth function of e G R defined by 

1 ix r 

n{e,x):= argr(- + y) + J dyuj{y,x), 
where the ambiguity in argr(i + y) is fixed arbitrarily and 



(6.34) 
(6.35) 

(6.36) 



d 1 ix 1 

..(y,x) :=- arg r(-y+ - + -) = - 



IX, 



, , 1 ix , , , 1 ..^ 



(6.37) 



It follows by means of ()5.10|) that uj{y,x) < 0, and one sees with the help of the asymp- 
totic expansion of and ()6.36|) (or directly from the asymptotic expansion of F) that 
lim^^-i-oo ^(e, a;) = ^oo. Therefore Q{e,x) decreases monotonically from +oo to —oo as e 
runs through the real axis. This implies that the set of solutions of ()6.34j) is bounded neither 
from below nor from above, and there are finitely many solutions in any finite interval. 
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In the foregoing derivation we assumed that k{X) ^ and fi;(A) 7^ cxd, but the conclusion 
remains vahd for these special values, too, as one can confirm by inspection of the respective 
conditions W^[p£;, V^i]o+ = and ^^[pe, V^2]o+ = 0. 

Since for A < the energy spectrum is not bounded from below, in physical applications 
one has to exclude those connection matrices V for which the angular Hamiltonian 
possesses a negative eigenvalue. 

7 Explicitly solvable cases 

We illustrate our inequivalent quantizations of the = 3 Calogero model by considering a 
few special cases which can be solved explicitly. Recall that our inequivalent quantizations are 
specified by the parameter k(A) (to each A < 1) for the radial part ()6.6p and the parameters 
in V fj3.18p for the angular part. In this section A will always be positive, and for simplicity 
we adopt the choice fi;(A) = for any < A < 1. Then the energy eigenvalues and the radial 
solutions have the form given by ()(i.3j) and ()(j.4j) /or all \ > 0. For the angular part, we 
consider the four cases, U = — 12, U = I2, U = ai and U = — o"i. We shall see, in particular, 
that the case U = ai admits a smooth limit for z/ — 1 (i.e., (7 — >■ 0) in which the system 
defined by Hrei ()2.3|) becomes the harmonic oscillator in two dimensions. For V ^ a\ the 
system resulting in the (7^0 limit can be interpreted as the harmonic oscillator plus an 
extra singular potential in the angular sector, which is supported at the points of S ()3.5|1 and 
manifests itself in the boundary condition. 

The spectrum of the angular Hamiltonian for f/ = ±l2, ±(Ti and the energy spectrum for 
\J = o\ are summarized by Figure 5 at the end of the section. 

7.1 The case U = -I2 

We begin by the 'Dirichlet' case V = — 12, which is in fact the standard choice of boundary 
condition that has been used since the introduction of the model [1 . As U = — 12 is one of 
the separating cases discussed in Section 4.2, we just recall fl4.21|) - ()4.24|) for the eigenstates 
of the angular part. With the eigenvalues, the solutions are 

6 

k=l 
6 

V^{<P) = J2^1<M)^ p = 2n + z., (7.1) 

k=l 

for n = 0, 1, 2, . . .. The arbitrary coefficients C^, k = 1, 2, . . . , 6 show that all levels have 
multiplicity 6. 

These solutions 77^, are then combined with the solutions for the radial part Rm,x{r) = 
Pm,x{r), where p^.A is given in ()6.4|) . to form the eigenstates for the entire system governed 
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by Hrei- Since A = is determined by ra, those states (modulo the normahzation constant) 
may be presented as (m, n = 0, 1, 2, . . .) 

^i^{r, 0) = i?„,A(r) r7;^(0), E^, = 2c (2m + 1 + 3(2n + 1 + z/)) , 

vl/^Jr, 0) = i?^,,(r) r^f^(0), = 2c (2m + 1 + 3(2n + u)) , (7.2) 



where 2c = y | see fl(j.3j) . 

On account of the multiphcity, one can choose the eigenstate in any particular represen- 
tation of the exchange-S'a (or more generally the Dq) symmetry group. For example, if we 
choose = (— l)*^"^^ for k = 2, ... ,6, then the resultant state rj^ for the angular part 
becomes a (+1) eigenstate of all the reflections P„ and hence it is bosonic. On the other hand, 
if we choose = Ci, then the resultant state r]^ becomes a (—1) eigenstate and hence it is 
fermionic. In contrast, for 77^ the choice = {—1)^^^C]_ makes it fermionic and = C\. 
makes it bosonic. At this point we recall that the basic solutions Vi^f^{(f)), i = 1,2, defined in 
()4.H1 are periodic in with period | and are symmetric with respect to = ^. Thus, if we 
introduce the sign factors, 

cos 30 sin 30 tan 30 

I cos 301 |sm30| I tan 301 

we can express these bosonic and fermionic eigenstates concisely over 5*^ \ 5 in terms of the 
basic functions. For instance, the bosonic states may be presented as 

<(0) = c(0) t;i,^(0), r^f (0) = t;i,^(0), (7.4) 

for the levels /i given, respectively, by ()7.H) . It is easily confirmed that the sign factor c(0) 
attached in ()7.4j) takes care of the sign conventions required for t]^ ^ together with the choice 
of the coefficients needed to provide the bosonic states. Similarly, the fermionic states are 
given by 

r/;f (0) = t(0) v,M, r^f (0) = .(0) t;i,^(0). (7.5) 
For the eigenvalues /i in 1)7.11) . one may use the relation 

+ + -;rin^3^j =j.^Cr(|coB3^|) (7.6) 

to replace Vi^^ in ()7.4|) and ()7.5|) with the Gegenbauer polynomial Cj". These bosonic and 
fermionic eigenstates recover the original solutions obtained by Calogero for = 3 |T]. 

7.2 The case U = I2 

The 'Neumann' case U = I2 can be solved analogously to the preceding 'Dirichlet' case; 
eq. ()4.21|) with ()4.25|) yields the angular solutions 

6 

k=l 
6 

= J2^lviM)^ /i = |2n + (1 - (7.7) 

k=l 
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for n = 0,1,2,..., with arbitrary coefficients C^, k = 1,...,6. The eigenstates and the 
eigenvalues for the entire system thus read 

vl>Kr, 0) = i?„,,(r) r/;^(0), i^l = 2c (2m + 1 + 3(2n + ! + (!- u))) , 
^Znir, <P) = RmAr) = 2c (2m + 1 + 3|2n + (1 - z/)|) . (7.8) 

As in the previous case, the eigenstates can be made bosonic or fermionic by choosing the 
coefficients C^, appropriately. Concise forms for these states are also available as 

<(0) = c(0) V2,M, 77^(0) = V2,M. (7.9) 
for the bosonic states, and 

<(0) =t(0) ^2,^(0), r^f(0) = si<l>)v2,M, (7.10) 

for the fermionic states. As before, by using the relation ()7.(i|l with u substituted by 1 — z/, one 
may replace ^2,/^ with the corresponding Gegenbauer polynomial Cj~'^ in the final expression 
of the solutions. 



7.3 The case U = ai 

The case ?7 = cxi is distinguished in the sense that it leads to the free connection condition 
(i.e., both the wave function and its derivative are continuous) in the limit u —>■ 1 where the 
singularity of the potential disappears. This 'free' case is one of the non-separating cases 
in which the spectrum consists of eigenvalues of both type 1 and type 2 eigenstates. The 
eigenvalues are determined by the spectral condition ()4.45|) . which now simplifies to 

-^ = 5J(r), (7.11) 

COSTTI/ 

using that the parameters in ()3.18|) are ^ = and B = 1. 

For type 1 eigenstates for which 9ft(r) = ±1, it is immediate to find the solutions for 
positive fi. For 3f?(r) = 1, these are 

l^ = 2n + l + iy and = |2n + (1 - z/)|, n = 0,1,2,..., (7.12) 

which correspond to the A_ and the case defined in fj4.48p and ()4.49p (see also ()4.23|) . 
()4.27|) ). Similarly, for 5R(r) = —1, we obtain 

/i = 2n + z/ and = 2n + 1 + {1 - u), n = 0,1,2,..., (7.13) 

which correspond to the B_ and the case. In fact, these are the solutions mentioned in 
()4.52j) - ()4.55j) . which arise since our parameters ()3.18|1 satisfy e*'-"'^^-' = 1 and e**^""'^-* = — 1. 
Thus, the type 1 solutions for the angular part are 

<+(0) = -c(0) ai{iJ,)v2,M), = 2n + 1 + (1 - z/), 

V^' (0) = ^(0) a2(/i)^i,M(0), /i = 2ra + 1 + z/, 

v!!H<P) = -bMv2,^{^), /i = |2n+ (1 - z/)|, 

?7f-(0) = s(0) b2{ix)viA(P). fi = 2n + u. (7.14) 
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Combining these with the eigenstates for the radial part, we obtain the type 1 eigenfunc- 
tions and energy eigenvalues for the entire system (m, n = 0, 1, 2, . . .): 

2c (2m + l + 3|2n+(l-z/)|), 
2c (2m + 1 + 3{2n + 1 + (1 - z/))) , 
2c(2m + l + 3(2n + z/)), 
2c(2m + l + 3(2n+l + z/)), (7.15) 

where the superscripts on "^mn specify the Dq representation similarly to ()4.57|1 . We observe 
that for U = ai the type 1 eigenstates 7]^+, rj^+ are basically the bosonic eigenstates ()7.9p 
admitted under U = I2, while r]^~ , rj^~ are the fermionic eigenstates ()7.5|) admitted under 
f/ = -12- 

Next, we turn to type 2 eigenstates for which 3?(r) = ±|- Since the spectral condition 
()7.1H) is analogous to the type 1 case, if we use 

A. N 1 

A u/j := — arccos - cosvrz/ 
^ ' -K \2 

we obtain the solutions, 

/i = 2n + 1 + A(i/) and /x = 2n + (1 - A(z/)), n = 0,1,2,..., (7.17) 
for 3?(r) = |, and 

/i = 2n + A(z/) and /i = 2n + 1 + (1 - A(z/)), n = 0,1,2,..., (7.18) 
for 3?(r) = -|. Note that 

\ < A(z^) < I A(z/) < u (7.19) 

for V in the range ()2.9|1 . We also remark that no solution with < is allowed for ()7.1H1 for 
both type 1 and type 2 eigenstates. 

The eigenfunctions associated with the type 2 eigenvalues can be constructed by the 
procedure of Section 4.3. Namely, one first obtains the eigenvector ()4.65|) with the aid of the 
projection operator 7r^(/i) in ()4.66|) . Then, taking into account ()4.30|) and ()4.43|) . one forms 
the eigenfunction rj^{(j)) in ()4.9|) out of the functions V±,n{4>)- Using the spectral condition 
()7.1H) and y{fi) = — 4iai(/i)a2(/i), z{fi) = — 4i6i(/i)62(Ai) for TTr{fi), and choosing the overall 
scale factor of the eigenvector appropriately, one arrives at 

^M,.(0) = -#t^i,.(0) + ^2,,(0), (7.20) 

where 

g(/x) = ^^2-2^ Ti-u + -) T{-u + -) T{u + ^^) T{u - /i). (7.21) 

This is valid for < < ^ in sector 1, and extension to the remaining half of sector 1 can be 
done by expressing in ()7.20j) the fj,^ on < < | in terms of the functions ?7±^^(0) defined by 
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vl/++(r,0) = i?„,,(r) 7/^-^(0), E++ = 

^- + (r, 0) = RmAr) ^mn = 

vl/+;(r,0) = (r)r^f(0), = 

'^^nij,(i)) = Rm,x{r)r]^-{(p), E-- = 



(7.16) 



()4.5|) . ()4.6|) . and adopting the resulting formula on the whole sector 1. Subsequent extension 
to sector k can be made in terms of the rotated functions ?7^^(0) defined in ()4.7p and the 
eigenvector for sector k which is given by t^^^ times the eigenvector for sector 1 (see ()4.30p 
and ()4.43|) ). Note that to each /i we have two solutions for r7^,r(0) on account of 53(r) = 
implying that each level is indeed doubly degenerate. 

To display the above eigenstates and their eigenvalues more systematically, let us use the 
notation introduced in the paragraph above ()4.6Up for the states belonging to the different 
type 2 representations of Dq. Thus the states rjl^T arise for 3fJ(r) = ^ and ?7/t,| for 3?(r) = — i, 
that is, for eigenvalues ()7.17|1 and ()7.18|) . respectively. We observe that, like in the case of 
type 1 states, each of the sets {ry^r} and {i]Ii,t} can be classified into two distinct series 

(2) — 

according to the difference in the non-integral part of /i. We introduce the notation rjii^r for 
the eigenstates ()7.20j) with fi = 2n + 1 + A{u) and ril^P~ for those with fi = 2n + (1 — A(z/)), 

~(2)— —(2)+ 

and similarly r/^,! with fi = 2n + A(z/) and ?7/j,t with = 2n + 1 + (1 — A(z/)). Combining 
with the solutions for the radial part, and adopting similar notation to specify the entire 
eigenstates containing the type 2 angular states, we obtain 

Ei^H = 2c (2m + 1 + 3{2n + (1 - A(z/)))) , 
E^2H = 2c (2m + 1 + 3{2n + ! + (!- A{u))) , 
Ei2)- = 2c(2m + l + 3(2n + A(z/))), 
El^l = 2c (2m + 1 + 3{2n + 1 + A(z/))) . (7.22) 

Note that fi can be recovered from the energy as 3/i = I - 2m - 1. The energy eigenvalues 
in ()7.15j) and ()7.22|1 provide the complete spectrum of the = 3 Calogero model defined by 
the Hamiltonian Hj-ei ()2.3|) under the angular boundary condition U = cxi. We mention that, 
for any the ground (the lowest energy) state is given by the type 1 state \I/qq^ possessing 
the energy = 2c (1 + 3|1 - z/|). 

Now, let us consider the harmonic oscillator limit z/ — » 1. Here, the functions in ()4.H) 
reduce to 

1 

= -sin3/x0, V2,f,{(p) = cos3/i0, (7.23) 

for sector 1, and we have 

ai(/i) = ^sin (^^j, a2(/x) = cos (^^), 

hifi) = 3cos (^), 62(/i) = -3/isin (^). (7.24) 

These are either zero or proportional to (—1)" for /i in ()7.14j) with z/ = 1, and hence the type 
1 states are basically given by the trigonometric functions ()7.23j) . Notice that in the u = 1 
limit, except for n = 0, the states ?7^+ and ?7^- are degenerate with eigenvalue /i = 2n, and 
similarly 1]^+ and ?7^~ are degenerate with /i = 2n + 1. These two pairs of degenerate states 
also share the same eigenvalue among themselves for T with r = 1 and r = — 1, respectively. 
Thus, one may form their linear combination to obtain the simpler set of eigenstates e^^^'^'^ 
for integers /i = 1, 2, . . .. These states have r = 1 for even and r = — 1 for odd. 



*^t;(^,0) = i?^,A(r)r/(^^(0), 
^^t(r,0) = i?„.,A(r)r^(^)-(0), 
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To find the type 2 states in the limit, we observe that A(z/) | for u 1, and that the 
factor in ()7.21|) reduces to 

= (7.25) 
4 sinvr/i 

Hence, for the eigenvalues /i in ()7.17|) . ()7.18p . the solution ()7.20p becomes ?7^,r(</') = e^'^^''^, 

which is valid for all sectors, where the signs ± correspond to Q{t) = for rjn^r = vl^^r^ 
and fjj^}-~ and to Q=(r) = =f^ for r/^^^ = ijl^}-' and r^^^l"^. 

Consequently, if we introduce k := 3fi (which yield integers for all as J/ — >■ 1), for u = 1 
both the type 1 and type 2 eigenstates can be combined to be presented together as 

77,^(0) :=e±''=^ A; = 0,1,2,.... (7.26) 

In the V = 1 limit, the complete set of eigenfunctions and eigenvalues of Hrei (j2.3j] is therefore 
furnished by 

^^fc(^ <P) = Rm,xir) %^(0), E^, = 2c (2m + 1 + fc) , (7.27) 
for m, k = 0,1,2, .. ., where it understood that both signs give the same if k = 0. In view of 
A = /c^, we see that 

i?™,A(r)=rV5-'L^(cr2). (7.28) 

The states (I7.26|l are the 27r-periodic eigenstates of the operator M in ()2.6|1 for (? = without 
singularity at the points of S ()3.5j) . Correspondingly, the eigenfunctions and the eigenvalues 
()7.27j] recover precisely the ones known for the harmonic oscillator in 2-dimensions (see, e.g., 
|22] for comparison). This shows that under the 'free' boundary condition U = ai our system 
is smoothly connected to the harmonic oscillator in the limit u —>■ 1. This is not the case 
for the previous two cases, U = — 12 and U = I2. Indeed, these become such systems in the 
u —* 1 limit in which 'two thirds' of the levels of the harmonic oscillator are missing and each 
level has multiplicity 6 (instead of 2) in the angular sector. 



7.4 The case U = -o\ 



The case \J = —cti, where we have ^ = and i3 = — 1 in ()3.18|) . can be dealt with analogously 
to the 'free' case U = ai. Indeed, the spectral condition ()4.45|) now reads 

^ = ^{r), (7.29) 
cosnu 

and hence, as a whole, the spectrum remains the same as that of the 'free' case. The only 
difference is that, because of the opposite sign in ()7.29|1 on the right hand side compared to 
(jZHH), the eigenvalues associated with the solutions are interchanged. For type 1 eigenstates, 
the interchange amounts to A± and B± B^. Thus, the angular solutions become 

^^+(0) = c(0) a2(/i)t;i,^(0), /i = 2ra + 1 + z/, 

<~ (0) = -t{<l>) ai{fi)v2,f,{(f)), fi = 2n + l + {l-u), 

^J^(0) = ^2(/i)fi,/.(0), /i = 2n + z/, 

Vu' (0) = -s(0) h{fi)v2M)^ fi = \2n + {l-u)\. (7.30) 
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Figure 5: Left: The angular 'eigenvalue-parameter' /i for the four explicitly solvable cases 
discussed, with u = 4/5. Boxes stand for 'case A' states and circles for 'case B' ones; the + 
or — within them shows the sign of 9?(r), a dot is applied when the sign is undetermined. 
The symbols H= and = indicate type 2 states with positive and, respectively, negative 3?(r). 
The superscript numbers display the multiplicities. Right: the energy spectrum (conveniently 
shifted and rescaled) for the case U — ui. 
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Accordingly, the type 1 eigenfunctions for the entire system are obtained from ()7.15|) with 
the interchange of eigenstates and eigenvalues as shown in ()7.30|) . 

The type 2 eigenstates acquire a similar change as observed for type 1 states. Explicitly, 
the solutions for the spectral condition are given by ()7.17|) and ()7.18p with the interchange of 
the cases 3?(r) = ^ and 3?(r) = — |, i.e., the states r]^"^^ and r^^^^ are swapped. Consequently, 
the type 2 eigenfunctions of the entire system are obtained from the solutions for the case 
f/ = (Ti by the corresponding interchange of eigenstates and eigenvalues ^. 

Finally, we mention that ii U = — cxi, then the system does not tend to the 2-dimensional 
harmonic oscillator as z/ ^ 1 , even though the spectrum reduces to that of the harmonic oscil- 
lator in this limit. This can be seen, for instance, by looking at the ground state wave function, 
Ro,\{f) (0) with = |1 — which has parity —1 under P„ for any n in disagreement with 
the +1 parity of the oscillator ground state. 

8 Conclusion 

In this paper we explored the inequivalent quantizations of the three-particle Calogero model 
in the separation of variables approach under the assumption ()2.9|1 on the coupling constant. 
Upon requiring the Dq symmetry, we found that the model permits inequivalent quantiza- 
tions for the angular Hamiltonian M fl2.6|) which are specified by boundary conditions of the 
form ()3.14p parametrized by a matrix U G f/(2) satisfying ()3.15|) . We showed that the angu- 
lar boundary conditions fall into the qualitatively different 'separating' and 'non-separating' 
classes, and it is possible only in the separating case to set the admissible wave functions to 
zero in all but one of the six sectors corresponding to the different orderings of the particles. 
Another important distinction was uncovered between the boundary conditions admitting 
and the ones not admitting a negative eigenvalue of M^, since in the former case the energy 
is not bounded from below, which is in general not permissible in physical applications. The 
properties that has an eigenvalue A < or that it possesses only eigenvalues A > 
are stable generically (in the sense of Section 5.4) with respect to small perturbations of the 
parameters of the 'connection matrix' U f)3.18|) . Our description of the inequivalent quanti- 
zations of the radial Hamiltonian ()6.2|) for < A < 1 is consistent with and complements the 
previous analysis [TUj . 

We classified the eigenstates of the Hamiltonian according to the irreducible represen- 
tations of the Dq symmetry group, and described also the induced classification under the 
exchange- 5*3 subgroup of Dq. If necessary in some application, one can consistently truncate 
the Hilbert space to a sector containing only the states of a fixed symmetry type. Our con- 
struction provides new quantizations also for the so-obtained truncated sectors, containing 
for example the states of 'bosonic' or 'fermionic' character with respect to the permutations 
of the particles. 

Our case-study illustrates the fact that inequivalent quantizations have very different 
properties in general, and external theoretical or experimental input is needed to choose 
between such quantizations. One possible criterion for the choice may be the smoothness of 
the model in the limit where the singularity of the potential disappears. Our solution for 
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U = ai mentioned in Section 7 shows that there indeed exists a distinguished quantization 
that meets this criterion. 

Of course, the present work can only be regarded as a 'theoretical laboratory' since most 
applications of the Calogero model use arbitrarily large particle number. It would be very 
interesting to extend our construction to the particle case, which would require understand- 
ing the possible self-adjoint domains of the partial differential operator Hq in under the 
assumption ()1.3|1 . For example, we wonder if an analogue of the explicitly solvable 'free' 
case that we found for N = 3 exists for general N. It would be also interesting to better 
understand the inequivalent self-adjoint domains of the Hamiltonian without adopting the 
separation of variables approach, starting directly from the minimal operator corresponding 
to the formal expression 

Note added. We learned after submitting the paper that the spectra of the self-adjoint 
extensions of TY^.a f)6.2|) have also been studied, for ^ < A < 1, in The method used in 
|2S1 is similar to that used in ^U], and the results are consistent with our results derived in 
Section 6 relying on a different, but equivalent, method. We thank P.A.G. Pisani for drawing 
our attention to this article. 
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A Remarks on the angular and radial Hamiltonians 

The characterization of the self-adjoint domains for the angular Hamiltonian described in 
Section 3 can be viewed as an application of the general theory of self-adjoint differential 
operators [71 El CHI- Nevertheless, it may be useful to present an elementary argument proving 
that the conditions in ()3.14p provide self-adjoint domains for M fl2.6|) . In this appendix we 
also wish to quote some theorems from [2] that imply the discreteness of the spectrum of the 
radial Hamiltonian Ti,. a ()6.2j) on any self-adjoint domain. 

It was mentioned in Section 3 that the self-adjoint domains for M arise as restrictions 
of the maximal domain T>i. Here, our aim is to show that the restriction T> G T>i defined 
by the conditions in ()3.14j) yields a self-adjoint domain, i.e., the restriction of M-p^ to "D is a 
self-adjoint operator. For this, it proves advantageous to rewrite ()3.14|) in the equivalent form 

UgBi+\'ilj) = bI,~\^) {UeeUi2),eeS) with bI;^\^) := Bei^) ±iB'sW . (A.l) 

Note that all the twelve 'boundary vectors' b\^\iI)) take independently all the possible 
vector values as ifj runs over Vi. One can see this by considering the boundary vectors 
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associated with the functions 

[c,v{ + c^v%, (A.2) 

where Ci, C2 G C, 6^ G 5 fl3.5|) is one of the singular points, and rj G C°°(S'^ \ S) is a function 
taking the constant value 1 on one side of 6' on a small closed interval and being zero on the 
other side of 9 as well as on both sides of the five other singular points. 

Next, let us point out that, for 'iip,rj G Vi, 



ees 

= J2 ((^^ B'oiv)) - (Boi^), B'M)) (A.3) 
ees 



ees 

where (■, ■) is the scalar product in L^(S'^) and (■, ■) is the scalar product in C^. Formula 
()A.3jl can be derived by partial integration using the identity 

which is valid on the domain of definition of the reference modes y?^ as a result of [y?^, yjg] = 1 
(see Section 3). It follows from (jA.lj) that, for tpji] &T>, the expression given by ()A.3j) vanishes 
(in fact, each term in the sum vanishes separately). This means that D is a symmetric domain 
within T>i, i.e., Mjy is a symmetric operator. To demonstrate that this domain is a self-adjoint 
one, it is enough to show that the vanishing of ()A.3|) for all 7/; G 7? with a fixed rj eVi implies 
that r]eV. 

We now choose two functions il'i,^2 G V such that, for a given 9 e S, 5j+^(V^i) and 
B^f^\ip2) form an orthonormal basis in [by (jA.lll . i?^ and Bg \ip2) then also form 

an orthonormal basis] and the boundary vectors at the other singular points are zero. If 

(V^,,M^,r/) - (M^,V^,,r/) = i[(i?J+)(^,), i?(+)(7;)) - {B\^\ij,), B'f\i^))] (A.5) 
is zero for A; = 1, 2, then we can write 

UeB^;\^) = Uo [(5j+)(^i),5j+)(r^))i?J+)(^i) + (<)(^2),fii+Hr^))fii+HV^2 

= {Bl;\^,),Bl;\r^))UeBl;\^,) + (5W(^2),i?r^(r]))t/.i?r^(^2) 

= {Bt\^'Pi),Bi-\^))Bi-\'^,) + {Bi-\,p,),Bl-\r^))Bl-\,p,) 

= B^-\^). (A.6) 

This implies that 77 G as required. 

The self-adjoint domains for the formal radial Hamiltonian Tir^x ()6.2j) can be treated simi- 
larly to the above, and this case is actually much simpler. Since the boundary condition ()6.6|) 
appears in several references [HllinillZI, we need not dwell on this point. We below summarize 
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the general results that imply the discreteness of the spectrum of the radial Hamiltonian on 
any of these self-adjoint domains. 

Recall that the 'discrete spectrum' of a self-adjoint operator consists of the isolated points 
of the spectrum that are eigenvalues of finite multiplicity, and the rest of the spectrum is 
called the 'essential spectrum'. (Note that the isolated points of the spectrum are always 
eigenvalues, and for ordinary differential operators all eigenvalues have finite multiplicity.) In 
the case of self-adjoint ordinary differential operators the essential spectrum is the same for all 
self-adjoint extensions of the 'minimal operator', and thus it can be assigned unambiguously 
to the underlying formal differential operator (see e.g. XIII. 6. 4 in pj). According to the 
statement of XIII. 7.4 0, the essential spectrum of the formal differential operator Tir^x on 
the interval (0, oo) decomposes as the union of the essential spectra of the operators of the 
same form on (0,xo] and on [xq, oo) for any xq > 0. The essential spectrum assigned to the 
interval [xq, oo) is empty by XIII. 7. 16 jZj, since the potential term of Hr,x tends to +oo as 
r — > oo. If A > |, then the potential also tends to +oo as r — 0, and the essential spectrum 
associated to (0, Xq] is therefore empty by XIII. 7. 17 7J. If A < 1, the same conclusion follows 
from XIII. 6. 12 in by using that the deficiency indices of the minimal operator on (0,a;o] 
are (2,2). By combining these, we see that the essential spectrum of the formal differential 
operator TCr,\ on (0, oo) is empty, and hence all of its self-adjoint versions have pure discrete 
spectra. 

The above arguments can be used to prove the discreteness of the spectrum of the radial 
Hamiltonian for any particle number A^, since H^^x given by ()1.5|) . ()1.7|) leads to the equivalent 
operator in L^(R_|_,(ir) (see eq. ()6.2|) ) 

N-2 2-N f)? (P N o o fl^ , ^ , ^ , o , o , . , 

nr,x ■■= r— o Hr,x o r— = - — -r^ + -mu\^ + — A^ - 2 A^ - 4 /r^ + \ r\ A.7 

2m ar^ 4 8m 

which must be a self-adjoint operator in (R_|_, c/r). For A^ = 3 the discreteness of the 
spectrum of the angular Hamiltonian also follows by similar reasoning, but for A^ > 3 in 
p.fi|) becomes a partial differential operator that would require a different treatment. 

B Representations of the symmetry group Dq 

The dihedral group Z^g admits four different 1-dimensional representations and two inequiv- 
alent 2-dimensional irreducible representations. This follows since the 12 elements of fall 
into 6 conjugacy classes as described in Figure IHl (with the notations in eqs. ()3.1|) - ()3.3|) ). and 
12 = 1 + 1 + 1 + 1 + 2^ + 22. 

The 1-dimensional (or 'type 1') representation of character with QiP E {+} is defined 
by assigning the parities g and p to the reflections and Pk {k = 1,2,3), respectively. 
Since the Pk generate the exchange- 5*3 subgroup of Dq, the representations with p = + can 
be called 'bosonic' and those with p = — can be called 'fermionic'. The character of the 
2-dimensional defining representation of Dq is denoted by x^"^^- The other 2-dimensional (or 
'type 2') representation is the tensor product of the defining representation and one of the 
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Figure 6: Character table of the group Dq. 



type 1 representations with character x"*" or x The type 2 representations of Dq remain 
irreducible (and become equivalent) when restricted to the 5*3 subgroups. 

For reference in the main text, note that the eigenvalues of T^^/g in the defining represen- 

27ri 

tation are —j and —J and in the other type 2 representation are j and J, with j = e~ . Indeed, 
this is a consequence of the relations X^^H^tt/s) = ^ = —J ~ J and X^'^^TZtt/s) = = J + J 
taking into account that the eigenvalues of TZ^^j,, must be sixth roots of unity. 



C The monotonicity of the function F 



A 



We here demonstrate that the function defined in fl5.2|) is strictly monotonically decreasing 
for n e with any m > in ()5.14p . as well as for /i G [0,/i[f ). 

Consider the logarithmic derivative of F4, 

'2-z/ + /i' 



U — jJ, 



(C.l) 



Remember that 

7/>(l - 2) = + TTCOtTTZ, (C.2) 

where ip{z) is strictly monotonically increasing on the positive real semi-axis, cotvr^ is de- 
creasing between two consecutive singularities. We can rewrite (jdljl as 



+ 
+ 



2 



1 - Z/ + /i 



vr cot — (z/ + yu) — vr cot — (1 — z/ + /i) 



(C.3) 



Using that | < z/ < | ()2.9|) . one sees that if 



/i > (2 - z/) 



(C.4) 
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then the arguments of the four ip functions in ()C.3|) as well as the contributions of the first 
two lines of this formula are positive. 

Referring to ()5.14|) . ()5.15|) for the notations, suppose now that 

/^m < < /^m, (C.5) 

where the function F4 is negative. By using the periodicity of cot, we find that 
cot — (z/ + yu) — cot — (1 — z/ + |u) = cot — (z/ + yU — 2m — 2) — cot — {1 — p + fi — 2m) > 0, (C.6) 

since for this range of /i 

< (z/ + /i-2m-2) < (l-z/ + /i-2m) < 2, (C.7) 

thanks to < | ()2.9j) . This proves that F'j^{fi) < 0, whenever Fa{ij) < (with yU > 0). 
For the 'positive branch' F^ifi) > 0, supposing that 

/i^ < /i < /i^+i, (C8) 

we obtain 

cot — (z/ + yu) — cot — (1 — z/ + /i) = cot 7r7 — cot vr-i? (C.9) 

with 

7:=^-(m + l), ^:=i^^-(m + l). (C.IO) 
The difference (l(y.9jl is negative, since 

0<^9<7<1, (C.ll) 

thanks to ()C.8|) and > | ()2.9|) . In order to combine the terms of different signs in ()C.3|) . we 
may use the following standard integral formulae: 

Hp)-Hq) = j^ — — dx, p,g>o, (C.12) 

7rcot7rg= / dX, < q < 1. (C.13) 

Jo 1 - X 

For /i in ()C.8|) . putting these into ()C.3|) using ()C.9|) . ()C.11|) and that now (1 — z/ + /x) > 0, we 

.^^<'')- ^"^WaY (C.14) 



FaM Jo 

with 



P(X) = (1 - X'"+^)(X^ - X'^)X-i + (1 - X^+'^+'"+^)(X-'^ - X-^), (C.15) 

as one can verify straightforwardly. Because of (jdllll . -P(X) < for < X < 1. This proves 
that < for /i in 
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The only case left to consider is 

0</i</i[J = (2-z/). (C.16) 

In this case the arguments of all four ip functions in (jdlj) are positive, and we may proceed 
with the aid of an integral formula relying on ()(I12j) . Similar to ()(I14jl . we now obtain 

P{X) = [X^" - X''^i)X^^ <0, (C.17) 
which completes the proof of our claim concerning the strictly decreasing nature of Fa- 
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